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Abstract 

We perform numerical studies of the running coupling constant an(p 2 ) and 
of the gluon and ghost propagators for pure SU(2) lattice gauge theory in the 
minimal Landau gauge. Different definitions of the gauge fields and different 
gauge-fixing procedures are used respectively for gaining better control over the 
approach to the continuum limit and for a better understanding of Gribov-copy 
effects. We find that the ghost-ghost-gluon-vertex renormalization constant is 
finite in the continuum limit, confirming earlier results by all-order perturbation 
theory. In the low momentum regime, the gluon form factor is suppressed while 
the ghost form factor is divergent. Correspondingly, the ghost propagator di- 
verges faster than 1/p 2 and the gluon propagator appears to be finite. Precision 
data for the running coupling o_r(p 2 ) are obtained. These data are consistent 
with an IR fixed point given by lim p _»o ^(p 2 ) = 5(1). 



PACS: 11.15.Ha, 12.38.Aw, 12.38.Lg, 14.70.Dj 

keywords: gluon propagator, ghost propagator, running coupling constant, 
SU{2) lattice gauge theory 



1 



1 Introduction 



The non-perturbative study of non-Abelian gauge theories is of great importance for the 
determination of infrared (IR) properties such as color confinement and hadronization. 
These properties are encoded in the low momentum behavior of Yang-Mills Green's 
functions. Derived from these Green functions, the so-called (renormalized) running 
coupling constant "a^Qo 2 )" plays an important role for phenomenological studies and 
model building. While at large momentum the running coupling decreases logarithmi- 
cally with momentum, it rapidly rises at the hadronic energy scale of several hundred 
MeVs thus signaling the breakdown of the perturbative approach. Nonperturbative 
studies of a R (p 2 ) may be carried out analytically using the Dyson-Schwinger equations 
(DSEs) and numerically through lattice simulations. 

Whereas the high momentum behavior of the running coupling is uniquely deter- 
mined and provided by perturbation theory, several definitions of the running coupling 
in the low momentum regime are possible. All of them match with the perturbative 
result at high energies. For example, the corrections to the Coulomb law of the static- 
quark potential may be used to define the running coupling PJ E] . Alternatively, the 
finite-size scaling has its imprint on the running coupling and may be used for a high- 
precision measurement of the coupling jSlHj- The approach adopted in [5] [7] is based 
on extracting the running coupling directly from a vertex function. For a recent review 
of lattice calculations for an(p 2 ) see jH]. 

In order to obtain Green's functions for the fundamental degrees of freedom, gluons 
and quarks, gauge fixing is necessary. Despite being gauge dependent, these Green 
functions play an important role for the phenomenological approach to hadron physics. 
Landau gauge is a convenient choice for gauge fixing for several reasons: First of all, it 
is a Lorenz-covariant gauge implying that 2-point functions only depend on the square 
of the momentum transfer. Secondly, the renormalization procedure is simplified since 
the ghost-ghost-gluon- vertex renormalization constant Z\ is finite, at least to all orders 
of perturbation theory. This result — obtained by Taylor jH] — is a particular feature 
of Landau gauge and allows another definition of the running coupling constant, which 
only requires the calculation of 2-point functions: Let Fr(p 2 , /i 2 ) and Jr(p 2 , /i 2 ) denote 
the form factors (for a renormalization point fi) of the gluon and the ghost propagator, 
respectively; the running coupling is then defined by 

a R (p 2 ) = a R (fi 2 ) F R (p 2 ^ 2 ) 4(p 2 ,/i 2 ) . (1) 

(see Subsection 12.71 below) . 

Due to its usefulness for the description of the physics of hadrons, the non-perturba- 
tive approach to low-energy Yang-Mills theory by means of the DSEs has attracted 
much interest over the last decade [TTH ITT] . The coupled set of continuum DSEs for the 
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renormalized gluon and ghost propagators in Landau gauge has been recently studied 
by several groups m an cases it was found that the gluon and ghost form 

factors satisfy simple scaling laws in the IR momentum range p <C 1 GeV 

F fl (pV) cx[p 2 ] a , J«(pV) oc [p 2 ] P (2) 

where the remarkable sum rule holds for the IR exponents a and (3: 

a + 2(3 = 0. (3) 

It is interesting that this result is rather independent of the truncation scheme under 
consideration. These exponents may be determined from lattice simulations, assuming 
the parameterization a = 2k and (3 = —k. Note that for k > this implies a divergent 
ghost form factor Jr(p 2 , /i 2 ) in the IR limit and a vanishing gluon form factor Fr(p 2 , /i 2 ) 
in the same limit. Also, since the gluon propagator is given by D(p 2 ) = F(p 2 )/p 2 , one 
gets that -D(O) is infinite or finite respectively if k < 0.5 or k > 0.5. In the second case 
one has -D(O) = for k > 0.5 and D(0) finite and non-zero for k = 0.5. The IR sum 
rule © also implies that the running coupling [defined in eq. (JoTj) below] develops a 
fixed point in the IR limit 

lima#(p 2 ) = a c = constant. (4) 

Note that this result is independent of the value of k as long as the IR sum rule (JBJ) is 
satisfied. 

The precise value of k as well as the fixed-point value a c depend strongly on the 
truncation of the Dyson-Schwinger tower of equations. In fact, depending on the 
truncation, one finds 0.3 < k < 1 in the four-dimensional case These studies 

vary in their vertex Ansatze, angular approximations of the momentum loop integral, 
and on the tensor structure considered. In Ref. ^2] a new class of truncation schemes 
has been introduced, which manifestly ensures the multiplicative renormalizability of 
the propagator solutions. In this truncation, although the exact values of k and a c 
depend on the details of the truncation of the DSE tower, the value of a c is constrained 
to 

2n/N c <a c < 8tt/N c (5) 

for SU(iV c ). 

An IR-finite gluon propagator and an IR-divergent ghost form factor j^Hl 

l3"T] are also obtained using numerical simulations in the minimal Landau gauge. The 
present numerical data for the ghost propagator indicate a value of k = —(3 smaller 
than 0.5, while for the gluon propagator it is still under debate if k = a/2 is equal 
to or larger than 0.5. In both cases large finite-size effects in the IR region make an 
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exact determination of these exponents difficult. This is particularly evident in the 
gluon propagator case [23 I2H I2H EU EH] , where one needs to go to very large lattices 
in order to have control over the infinite-volume extrapolation. Our present data are 
consistent with k of the order of 0.5 (SUESj- 

Let us stress that in the minimal Landau gauge, which is the gauge-fixing condition 
used in numerical simulations (see Section 12731 below) . the gauge-fixed configurations 
belong to the region of transverse configurations, for which the Faddeev- Popov operator 
is non-negative. This implies a rigorous inequality [36J-[HB] for the Fourier components 
of the gluon field and a strong suppression of the (unrenormalized) gluon propagator 
in the IR limit. At the same time, the Euclidean probability gets concentrated near 
the border of this region, the so-called first Gribov horizon, implying the enhancement 
of the ghost propagator at small momenta [SB] • A similar result was also obtained by 
Gribov in [HHj . 

Taylor's finding is based upon the Faddeev-Popov quantization, thereby ignor- 
ing the effect of Gribov copies, which are certainly present in an intrinsically non- 
perturbative approach. The goal of the present paper is to confirm Taylor's result (Z\ 
is finite to all orders of perturbation theory) by the non-perturbative approach provided 
by lattice simulations. In addition, a thorough study of the gluon and the ghost form 
factors is performed. A focal point is the IR limit of the running coupling constant. 
Support for the existence of the fixed point is found, and a first estimate of a c is pro- 
vided from extensive lattice simulations. We present two sets of simulations, carried 
out respectively in Sao Carlos and in Tubingen, employing different definitions of the 
gauge fields and different gauge-fixing procedures. We believe that the comparison of 
these two formulations strengthens the importance of our findings. 

The paper is organized as follows. In Section El we describe the lattice approach to 
Yang-Mills Green's functions. Section |3] contains the numerical setup for the simula- 
tions carried out in Sao Carlos and in Tubingen. In Section 0] we report our data for 
the gluon and ghost propagators and for the running coupling constant. Conclusions 
are left to the final section. 

Preliminary results have been presented in [33J EHl IS] ■ 

2 The lattice approach to Green's functions 

In this section we explain the two lattice setups used for the numerical evaluation of 
the gluon and ghost propagators. We also recall the definition of the running coupling 
constant considered in Ref. [12], which can been evaluated using these propagators. 
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2.1 Gluon field on the lattice 

The action S of the continuum SU (2) Yang-Mills theory is formulated in terms of the 
field strength 

F; v [A]{x) = d,A a v {x) - d u A;(x) + g e abc A b ^x) Al{x) (6) 
and is given by 

S = \J d^x F^[A](x) F^[A](x) . (7) 

Here go is the bare coupling constant and A a ^(x) is the continuum gauge field. 

On the lattice the dynamical fields are SU(2) matrices U^x), which are associated 
with the links of the lattice, and the Wilson action is given by 

S = & E 1 " \^cP,u(x), (8) 

X,fJ,>U 

where the plaquette is defined as 

P^(x) = U^x) U v {x + e„) Ul(x + e u ) Ut{x) (9) 

and is a unit vector in the positive /i direction. (Note that the trace extends over 
color indices only.) The Wilson action is invariant under the gauge transformation 

U$(x) = Q{x) U^x) Q\x + e M ) , (10) 

where Q(x) are SU(2) matrices. The link variables U^x) may be expressed in terms 
of the continuum gauge field A^x) by making use of the relation 

U^x) = exp [iag Al(x)t b ] , (11) 

where a is the lattice spacing, t b = <r b /2 are the generators of the SU(2) algebra and 
a b are the Pauli matrices. One can check that in the naive continuum limit a — > the 
Wilson action reproduces the continuum action in eq. (jZJ if 

13 = A/g 2 = l/(7ra ) , (12) 
where ao is the bare coupling constant (squared). 

For the gauge group SU(2), the link variables U^(x) can be given in terms of (real) 
four-vectors of unit length 

U^x) = uJ(x)]L + zr? M (z)a , [ul(x)] 2 + [u^xf = 1 , (13) 
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where 1L is a 2 x 2 identity matrix. By defining the lattice gluon field A b (x) as 



one obtains 

A%x) = 2ul(x) = ag Al(x) + 0(a 3 ) (15) 

in the naive continuum limit a —>■ 0. This definition has been used for the numerical 
simulations done in Sao Carlos. 

Note that the gluon field A b (x) defined above changes sign under a non-trivial center 
transformation Z2 of the link fields U^x) —>■ —U^x). Recently, another identification 
of the gluonic degrees of freedom in the lattice formulation was proposed [2E|- In this 
case, one first notices that the gluon field in continuum Yang-Mills theories transforms 
under the adjoint representation of the SU (2) color group, i.e. 

,aed 

A a ;(x) = O a \x)Al{x) + — O ec (x)d,O dc (16) 
O a \x) = 2tr c [n{x)t a tf{x)t b ] , (17) 

where Q(x) G SU(2) is a gauge transformation of the fundamental quark field and 
O ab (x) G 50(3). In view of the transformation properties in eq. tj!6j). one can iden- 
tify the continuum gauge fields A a (x) with the algebra- valued fields of the adjoint 
representation 

Uf{x) = {e W [ag Al{x)i b ]} cd , (18) 

where i b ac = e abc and the total anti-symmetric tensor e abc is the generator of the £77(2) 
group in the adjoint representation. On the lattice, the adjoint links U^ b (x) are obtained 
from 

Uf{x) = 2tr c [U,(x)t c Ul(x)t d ] (19) 
and the gluon field A®(x) is given by 

A%x) = 2u%x)ul{x), (20) 

without summation over /i on the right-hand side. By expanding eq. ()18j) in powers of 
the lattice spacing a and by using eqs. (fTB^) and (JTHJ) one obtains 

A%x) = ag,A%x) + 0(a 3 ) . (21) 

Clearly, the representation ()20j) is invariant under a non-trivial center transformation 
Ufx{x) — > — U/j,(x). This discretization of the gluon field has been used for the simula- 
tions done in Tubingen. 
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It is well known that different discretizations of the gluon field lead to gluon prop- 
agators equivalent up to a trivial (multiplicative) renormalization [12] IIS]- Also, this 
proportionality constant between different discretizations of the gluon propagator may 
be (partially) explained as a tadpole renormalization jUJIH]. We point out, however, 
that it is useful to disentangle the information carried by center elements and coset 
fields, defined above, when the vacuum energy is investigated. In particular, it was 
found that — in the continuum limit — the center elements provide a contribution to 
the gluon condensate [IE1 EHj ■ 



2.2 Tadpole improved gluon fields 

The relation between lattice and continuum gluon fields relies on the expansion [see 



1 + iag A b (x)t b + ... 



(22) 



where the ellipses denote higher order terms in the bare coupling constant go- The 
artificial contributions of these higher order terms to loop integrals are called "tadpole" 
terms [SUj. These terms are only suppressed by powers of g\ and are generically large 
in simulations using moderate (3 values. 

In order to remove the tadpole contributions from the observable of interest one 
can redefine the relation between the link matrices and the continuum gluon field by 
using 



UQ,L 



1 + iag A b (x)t b + 



where Mo,l is given by the "meanfield" value 

'tr. 



2 U,(x 



(for arbitrary fx) 



(23) 



(24) 



with the links U^(x) fixed to the Landau gauge. Equivalently, one can use (HO] a gauge 
invariant definition of the tadpole factor given by 



Uo,p 



^ 2 ^ D P' V ^ X ' 



1/4 



(25) 



where P^ u (x) is the plaquette defined in eq. Q. Thus, the use of tadpole improvement 
in the case of the standard definition of the lattice gluon field ([THJl gives 



Afa) = agoAfa) + 0(a 3 ) = 2u%x) / 



(26) 



In the case of the coset definition (j2~Uj) of the gauge fields, the tadpole factors are 
expressed in terms of the expectation values of the adjoint link and of the adjoint 
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plaquette and are given by 




Thus, the tadpole improved relation between the continuum gauge field A (x) and the 
link matrices is given in this case by 

A%x) = ag A b ^x) + 0(a 3 ) = 2 / < P . (29) 

In the sections below we will stress the effect of tadpole improvement on the various 
quantities considered in this work. 

2.3 Minimal Landau Gauge 

The gluon and ghost propagators depend on the choice of the gauge. In order to 
maintain contact with the Dyson- Schwinger approach and the results presented in the 
Introduction, we consider the so-called minimal (lattice) Landau gauge. This gauge 
condition is imposed by minimizing the functional 

S fix [fi] = -£tr c t#(x), (30) 

where Uff(x) is the gauge-transformed link (jlOj) . This minimizing condition corresponds 
to imposing the transversality condition 

(A ■ A) b (x) = A % x ) ~ A % x ~ e v) = V 6 and x , (31) 

which is the lattice formulation of the usual Landau gauge-fixing condition in the 
continuum. Let us notice that the condition (f3Tj) is exactly satisfied by the lattice gauge 
field only if the standard discretization (JT5j) is considered, while for the discretization 
given in (j2*U|) the above result is valid up to discretization errors of order a 2 . However, in 
both cases, the gauge-fixing condition (J3Uj) implies that the continuum Landau gauge- 
fixing condition d- A = is satisfied up to discretization errors of order a 2 . In practice, 
we stop the gauge fixing when the average value of [(A • A) b (x)} 2 is smaller than 10 -12 . 

The minimizing condition (j3*Uj) also implies that the Faddeev-Popov matrix is pos- 
itive semi-definite. In particular, the space of gauge-fixed configurations {Uf?(x)} lies 
within the first Gribov horizon, where the smallest (non-trivial) eigenvalue of the 
Faddeev-Popov operator is zero. It is well known that, in general, for a given lat- 
tice configuration {U^x)}, there are many possible gauge transformations Q(x) that 
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correspond to different local minima of the functional ()3U|) . i.e. there are Gribov copies 
inside the first Gribov horizon [SniETj. Thus, the minimizing condition given in eq. 
(}30|) is not sufficient to find a unique representative on each gauge orbit. A possible so- 
lution to this problem is to restrict the configuration space of gauge- fixed fields Uj}(x) 
to the so-called fundamental modular region i-e. to consider for each configuration 
{U^x)} the absolute minimum of the functional (JHUj) . From the numerical point of 
view this is a highly non-trivial task, corresponding to finding the ground state of a 
spin-glass model 51. On the other hand, if local minima are considered, one faces the 
problem that different numerical gauge-fixing algorithms yield different sets of local 
minima, i.e. they sample different configurations from the region delimited by the first 
Gribov horizon (see for example jH] and references therein). This implies that nu- 
merical results using gauge fixing could depend on the gauge-fixing algorithm, making 
their interpretation conceptually difficult. 

In the simulations done in Sao Carlos, a stochastic overrelaxation algorithm |52j-[E%j 
was used. The simulations done in Tubingen have employed the simulated annealing 
technique described in detail in Ref. [20] • The problem of Gribov copies was not 
considered in either case. Even though neither method is able to locate the global 
minimum of the gauge-fixing functional, i.e. to restrict the gauge-fixed configuration 
space to the fundamental modular region, a comparison of the propagators obtained 
using the two methods can provide an estimate of the bias ( Gribov noise) introduced by 
the gauge-fixing procedure. From this comparison we have found that the data for the 
propagators are rather insensitive to the particular choice of gauge-fixing algorithm, 
suggesting that the influence of Gribov copies on the two propagators (if present) is at 
most of the order of magnitude of the numerical accuracy. For the gluon propagator 
this result is in agreement with previous studies in Landau gauge for the SU(2) and 
SU(3) groups in three and four dimensions [21 EH EE]- A similar result has also 
been obtained for the gluon propagator in Coulomb gauge jSHl- On the other hand, 
a previous study of the ghost propagator in SU(2) Landau gauge [3T| has shown a 
clear bias related to Gribov copies in the strong-coupling regime. In particular, data 
(in the IR region) obtained considering only absolute minima have been found to be 
systematically smaller than data obtained using local minima. This result — which 
has been recently confirmed in jSTj ~ can be qualitatively explained. In fact, as said 
above, the smallest non-trivial eigenvalue X min of the Faddeev-Popov operator goes to 
zero as the first Gribov horizon is approached. At the same time, one expects that 
global minima (i.e. configurations belonging to the fundamental modular region) be 
"farther" away from the first Gribov horizon than local minima. Thus, the absolute 
minimum configuration should correspond to a value of Xmin larger — on average - 
than the value obtained in a generic relative minimum. 1 Since the ghost propagator is 
given by the inverse of the Faddeev-Popov matrix (see Section EUD below), this would 

1 This was checked numerically in Ref. |58|. 



9 



imply a smaller ghost propagator (on average) at the absolute minimum, as observed in 
Refs. [HU GUj • The analysis carried out in these references has shown that Gribov-copy 
effects are visible only for the smallest nonzero momentum on the lattice, at least for 
the lattice volumes considered, which are still relatively small. As explained below (see 
Section |Hn|) ) in our analysis we have not considered the data points corresponding to 
the smallest momenta. 



2.4 Gluon propagator 

The continuum gluon propagator in position space is given by 

K( X ~V) = (A a ,(x)Al(y)). (32) 
Correspondingly, one can consider the (position space) lattice gluon propagator 

V*{x-y) = (A;(x)Al(y)), (33) 

where A a Ax) is one of the lattice discretizations of the continuum gluon field discussed 
above [see eqs. ([To]) and (pjj)]. At the leading order a these two quantities are related 
by 

gl a 2 D*{x -y) = V%(x - y) / u\ P , (34) 

where u ,p is the tadpole factor given in eq. (|23j) [respectively eq. (J2~%j) ] when considering 
the lattice gluon field defined in eq. (|T5j) [respectively eq. (j2~Uj) ]. 

The lattice gluon propagator in momentum space is obtained by evaluating the 
Fourier transform 

vt(v) = ^XX^ - y) ex p I**- ( £ - *)] ' h = ^ ' (35) 



x,y >* 



where n M labels the Matsubara modes in the [i direction, is the number of lattice 
points in the same direction, x = xa, y = ya and V is the lattice volume. In order to 
minimize discretization effects , we consider the gluon propagator as a function of 
the lattice momentum p with components 

p, = 2 sin fe) . (36) 
It is also useful to introduce the lattice gluon form factor ^(p 2 ), defined as 

2>(p) = ^, m = \Y, v i» . ( 3 7) 

which is a measure of the deviation of the full propagator from the free one. Note that, 
in Landau gauge, the propagator is diagonal in color space and transversal in Lorentz 
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space. The transversality condition ()31|) implies [3T] that T>(0) is not given by T>(p) at 
p — 0. In fact, for p = the previous equation becomes ©(0) = 1/12^ 2^(0). 

In order to evaluate numerically the gluon propagator in momentum space it is 
useful to employ the formula [3~T] 



9V ^ 



2 ^M^) cos (P- 

2 ^(x) sin (p • x 



+ 

In fact, by expanding the previous equation we obtain 

V ^ = oi7 E E ( A » A W ) cos \p ■ ^ 



(38) 



(39) 



a,fi x,y 



which is directly related to eq. (JSEJ)- 

One can also evaluate the form factor F(p 2 ) directly |2Hj- To this end we can 
consider, without any loss of generality, a momentum transfer parallel to the time 
direction p = (0, 0, 0,p 4 ) and define 



(40) 



&tA^{x) = A^(x + e 4 ) 
where e 4 is the unit vector in the time direction. The form factor is then obtained from 

2 

1 / 

s2> 



HP 



-E 



a, 



+ 



^A^(x) cos(p-x) 

x 

J2 A tA a ^(x) sin(p-x) 



(41) 



By expanding the previous formula one can verify that the free part 1/p 2 is canceled 
exactly. This strongly suppresses the statistical noise in the high-momentum regime. 
Here, we will present results that directly address the gluon propagator and the 
gluon form factor (jSJ), evaluated respectively in Sao Carlos and in Tubingen. 



2.5 Ghost propagator 

The ghost propagator G ab {p) is uniquely defined once the gauge-fixing functional (|3T)|) 
is specified. In fact, if we write the gauge-fixing matrix as 

Q(x) = exp[i6 a {x)t a ] , (42) 
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with t a defined as in Section 12.11 then the gauge-fixing functional can be expanded 
with respect to the angles 9 a (x) and, at any local minimum of S^, we obtain 

s fi x = So + 0a ^ K eb (y) + °( d3 ) > ( 43 ) 

x,y a,b 

where is the so-called Faddeev-Popov operator. Note that the linear term in 9{x) is 
absent by virtue of the minimizing gauge-fixing condition [see eqs. (|30jl and l)31|) ]. The 
expression of the Faddeev-Popov operator in terms of the gauged-fixed link variables 
can be found in jHEJ eq. (B.18)]. Note that the matrix M^ b obtained in this way is a 
lattice discretization of the continuum Faddeev-Popov operator (—d + A) ■ d and that 
this discretization yields automatically the standard discretization A b (x) for the gluon 
field given in eq. ([T5]). 

The lattice ghost propagator Q ab {p) is provided by the inverse Faddeev-Popov op- 
erator M^y. Due to translation invariance, the lattice average of the inverse operator 
depends only on (x —y). Thus, in momentum space we have 

<3 a \P) = - E<( M-1 )0 ex P Hp • ( £ - v)\ ■ ( 44 ) 

Since the matrix depends linearly on the link variables U^(x), tadpole improvement 
applied to the ghost propagator implies a rescaling 

g ab ( P ) - g ab (p) u ,p . (45) 

Thus, at the leading order a one has 

a 2 G ab (p) = g ab (p)u 0:P , (46) 

where G ab (p) is the continuum ghost propagator in momentum space. 

The asymptotic behavior of the ghost propagator G ab (p) is known from perturba- 
tion theory: it decreases as 1/p 2 with additional logarithmic corrections. The 1/p 2 
behavior is inherited from the free-theory case. The non-trivial information on the 
ghost propagator is therefore encoded in the (continuum) form factor J(p 2 ), which is 
defined by 

i :. \ a"'' f ' i :. \ — a'"'' ' 

yielding the lattice form 



G ab {p) = 5 ab G(p) = 5 ab ^l , (47) 



g ab {p) = 5 ab ^p- . (4? 



p* 

Numerically, the lattice ghost propagator can be obtained by inverting the Faddeev- 

rab 



Popov matrix M£ b . In the numerical simulations carried out in Sao Carlos this has 
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been done using a conjugate-gradient algorithm. On the contrary, in the simulations in 
Tubingen the (lattice) ghost form factor J {ft 2 ) has been evaluated directly. To this end 
one can consider the following set of linear equations (for a given set of link variables 
U) 



K b y [U] u\y) = n a {cos [ft ■ (x - e,)} - cos [ft ■ x}} 



y,b 



Ky[U] v\y) = n a {sin [ft-(x- e,)} - sin [ft ■ x}} 

y,b 



(49) 
(50) 



where n a is an arbitrary unit vector which specifies the components of the ghost prop- 
agator under investigation. We are considering momenta with a non-zero component 
only in the /i direction. In fact, by solving these equations for u b (y) and v b (y) and 
using trigonometric identities we find that the ghost form factor is given by 



j (ft 2 ) = — ^2 ( { cos \p-(v- e v)\ - cos \P ■ y\ } nb u \y) 



(51) 



+ {sin [ft - (y- e M )] - sin [ft ■ y]} n b v\y) 



x,y a,b 



4 sin^ I ^ 



cos [ft ■ (x 



(52) 



Note that, with our choice of momenta, the lattice momentum squared is given by 
p 2 = 4 sin 2 (j3 M /2), implying that the free part 1/p 2 of the ghost propagator exactly 
cancels out and we are left with the form factor J (ft 2 ). The set of equations (|4lH) - (|5T)|) 
has been solved using a bi-conjugate gradient method for matrix inversion. 



2.6 Renormalization 

Renormalization of Yang-Mills theories in four dimensions implies that the bare cou- 
pling acquires a dependence on the ultraviolet (UV) cutoff kuv given by 

a -> a (A uv / A scale ) . (53) 

Thereby, the bare coupling constant is no longer the theory's parameter. The Yang- 
Mills scale parameter A sca i e takes over the role of the only parameter of the theory. 
In the context of (quenched) lattice gauge simulations the string tension a is widely 
used as the generic low-energy scale. In this case, the cutoff dependence of the bare 
coupling is implicitly given by the (3 dependence of a a 2 ((3) where (3 is related to the 
bare coupling in (fT^ and Auv = ^/ a (P)- In Section l3~3l we will derive this relation 
from lattice data (obtained in Ref. |60j). 
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In addition, wave-function renormalization constants develop a dependence on 
A uv / A sca i e . The lattice bare form factors of the previous subsections, Tb and J B , 
are related to their continuum analogues (for very large (5) by 

F B = Fb , 4 , Jb = Jb «o,p • (54) 
4m o,p 

These form factors depend on the momentum p 2 and on the UV cutoff Ajjy (given in 
units of the string tension) or, equivalently, on the lattice coupling f3 (see Section 14.31 
below). Thus, we can write 

T B = Tb(p 2 ,P) , J B = Jb(p 2 ,P) • (55) 
The renormalized form factors are obtained upon multiplicative renormalization 

F R (p 2 ,^ 2 ) = Zj 1 (p, f i)F B (p 2 ,f3) (56) 
J R (p 2 ,fi 2 ) = Z^((3,fi) J B (P 2 ,P) (57) 

using the renormalization conditions 

F R (fx 2 ,fi 2 ) = 1, Jp(/i 2 ,/i 2 ) = 1. (58) 

(Notice that tadpole renormalization does not affect the calculation of F R , J R but may 
be useful for the determination of the renormalization constants Z 3 , Z 3 .) 

Clearly, similar relations hold also for the bare and renormalized gluon and ghost 
propagators. In practice, the multiplicative renormalizability of the theory implies 
that a rescaling of the data for each (3 value (independently of the lattice momentum) 
is sufficient to let the form factors !Fb(p 2 ,P) and Jb(p 2 ,/3) - - or equivalently the 
corresponding propagators — fall on top of a single curve describing the momentum 
dependence of the corresponding renormalized quantity. 



2.7 Running coupling constant 

Of great importance for phenomenological purposes is the running coupling strength 
a n{p 2 ) considered in Ref. [T2j. In particular, this strength enters directly the quark 
DSE and can be interpreted as an effective interaction strength between quarks [61. 
This running coupling strength is a renormalization-group-invariant combination of the 
gluon and ghost form factors. In order to derive this combination we can start with 
the definition of the ghost-ghost-gluon- vertex renormalized coupling strength 

2 Z 3 (/3,fj) Z 2 (P,fi) 

a R {ji ) = ~ «o(A[/vO , (59) 

Zf(l3,fi) 
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where Zi(/3,/j,) is the ghost-ghost-gluon- vertex renormalization constant. In lattice 
simulations, the UV-cutoff is related to (3 by A uv = ir/a(f3), where a is the lattice 
spacing. Using eqs. (J56|) and (jBTjl we can express the renormalization constants Z 3 ((3, n) 
and Zs(j3,fi) in terms of the bare and renormalized form factors yielding 

a R (^) F R tf,f) = ■ ( 6 °) 

Note that the left-hand side of this relation is finite and independent of j3 by construc- 
tion and that the right-hand side depends on the renormalization scale /i only through 
the ghost-ghost-gluon- vertex renormalization constant Zi(j3,fi). It was found more 
than twenty years ago by Taylor [H] that (in the continuum) Z\ is finite, independent 
of /i, at least to all orders of perturbation theory. This finding will be confirmed by our 
lattice studies below. In this case, the right-hand side of eq. (JoTij) is thus independent 
of fi. Then, by choosing = a/p 2 and using the renormalization conditions (p)Hj) . we 
find the final expression for the running coupling strength, i.e., 

a R ( P 2 ) = a R (ji 2 ) F R (p 2 ,fi 2 ) J 2 R (P 2 ^ 2 ) • (61) 

Finally, let us notice that eqs. (J34"|) and (|47)|) imply that tadpole renormalization 
does not affect the renormalized coupling defined above. 



3 Details of the numerical simulations 
3.1 Setup 

All our simulations used the standard Wilson action for SU(2) lattice gauge theory 
in four dimensions with periodic boundary conditions. In order to check finite-volume 
effects and verify scaling we consider several values of (3 and of the lattice volumes 
V = x N t . The dependence of the lattice spacing on (3 can be inferred from a 
calculation of the string tension a in lattice units. Here, we will use the data for era 2 
reported in Ref. [60] and a linear interpolation of the logarithm of these data where 
needed. 2 A value of the lattice spacing in physical units was obtained using the value 
a = [440 Me V] 2 for the string tension. We used iV con f independent configurations for 
the numerical evaluation of the propagators. 

Computations in Sao Carlos were performed on the PC cluster at the IFSC-USP 
(the system has 16 nodes with 866 MHz Pentium III CPU and 256 MB RAM memory). 
All runs in Sao Carlos started with a random gauge configuration and for thermalization 
we use a hybrid overrelaxed (HOR) algorithm. The total computer time used for the 
runs was about 50 days on the full PC cluster. In Table Q we report the parameters 

2 For (3 = 2.15 an extrapolation of these data was necessary in order to obtain aa 2 . 
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(3 2.2 


2.3 


2.4 


2.5 2.6 


2.7 2.8 


aa 2 0.220(9) 


0.136(2) 


0.071(1) 


0.0363(3) 0.018(1) 


0.0103(2) 0.0055(3) 



Table 1: Simulation parameters of the runs in Sao Carlos. Data from Ref. 60J were 
used to obtain the lattice spacing in units of the string tension. Error bars (in paren- 
theses) come from propagation of errors and indicate one standard deviation on the 
last significant digit. The lattice volumes V and the number of configurations iV con f 
considered are discussed in the text. 



p 


2.15 


2.2 


2.3 


2.375 


2.45 


2.525 


aa 2 


0.280(13) 


0.220(9) 


0.136(2) 


0.083(2) 


0.0507(8) 


0.0307(5) 


iV 3 x N t 


16 3 x 32 


16 3 x 32 


16 3 x 32 


16 3 x 32 


16 3 x 32 


16 3 x 32 


IVconf 


200 


200 


200 


200 


200 


200 



Table 2: Simulation parameters of the runs in Tubingen. Data from Ref. (SU] were 
used to obtain the lattice spacing in units of the string tension. 



used for the simulations in Sao Carlos. For each (3 value three different lattice volumes 
were considered, i.e. V = 14 4 , 20 4 and 26 4 . For the lattice volume V = 14 4 (respectively 
V = 20 4 and 26 4 ) and for each (3 value we produced iV con f = 500 (respectively 150 and 
50) configurations. 

Computations in Tubingen were carried out at the local PC cluster where 4 — 12 
nodes (1 Ghz Athlon) were used. The simulation parameters for the runs in Tubingen 
are listed in Table El 

Table El lists the differences between the Sao Carlos and the Tubingen approach. 
Reconstructing the continuum gauge field from the link fields in different manners 
[compare eq. ()15|) with eq. ()20j)] provide insight into the discretization errors. Employ- 
ing different gauge-fixing algorithms points out the effect of the Gribov ambiguities on 
the propagators. 





Sao Carlos 


Tubingen 


definition of gauge fields 


fundamental representation 


adjoint representation 


gauge fixing 


iterative stoch. overrelaxation 


simulated annealing 


number of lattice points 


finite-size control 


fixed 



Table 3: Differences between the Sao Carlos and the Tubingen numerical approach. 
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3.2 Determination of renormalization constants 

In order to obtain the renormalized propagators and form factors, one needs to eval- 
uate the renormalization constants (/3,/i) and Z 3 " 1 ( ( 5,/x) defined in eqs. f!56j) and 
(|57|) . respectively. Multiplicative renormalizability implies that one can "collapse" data 
obtained at different (3 on a single curve. This can be done by using the matching tech- 
nique described in detail in Ref. jH2 Sec. V.B.2], For instance for the gluon form factor 
this is equivalent to considering the quantity 



where the factor Z% /x) for each (3 is obtained from the matching technique. The 
renormalization point, i.e. the \x dependence, comes into play when the "single" curve 
is rescaled to satisfy the condition 



The same procedure is applied to the ghost form factor. For our analysis we considered 
a renormalization scale of \i = 3 GeV. 

For the Sao Carlos data we checked for finite-size effects before applying the match- 
ing technique. 3 In particular, by comparing data at different lattice sizes and same (3 
value, we find (for each (3) a range of momenta for which the data are free from finite- 
volume corrections. We then perform the matching using data for these momenta and 
V = 26 4 . 

3.3 Asymptotic scaling 

As said in Section 12. ti\ the cutoff dependence of the bare coupling is implicitly given 
by aa 2 (/3). Thereby, the string tension a serves as the fundamental energy scale. 

In principle, perturbation theory predicts the (3 dependence of aa 2 for the regime 
a <C \j yfo. In practice, large deviations of the measured function aa 2 {(3) from the per- 
turbative scaling are observed in the scaling region, which corresponds to the interval 
(3 G [2.15, 2.8] in our case. Clearly, the relation between the "measured" values for aa 2 
and perturbation theory is highly important for a careful extrapolation of the lattice 
data to continuum physics. One goal of the present paper is to present this relation. 

For this purpose, we perform a large f3 expansion of the lattice spacing in units of 
the string tension, i.e., 



F R (p 2 ,^ 2 ) = Z^\(3,^)F B ((3,p 2 ), 



(62) 



F R (fi 2 ,fx 2 ) = 1 . 



(63) 




(64) 



3 For more details see (321 Sec. Ill 
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Figure 1: The string tension in units of the lattice spacing: lattice data from [60] and 
the fit using eq. (JMj) . 



The first two terms on the rhs of ()64j) are in accordance with 2-loop perturbation 
theory. The term d/f3 represents higher-order effects and the term c is a dimensionless 
scale factor to the string tension. Parameters c and d are determined by fitting the 
formula (j^U) to the lattice data reported in Ref. jHO]. Using only data for (3 > 2.3 we 
obtain 

c = 4.38(9) , d = 1.66(4) , ^jd.o.f. = 0.62 . (65) 

The corresponding fit is shown in Figure ^ It appears that the truncation of the series 
(}64|) at the 1//? level reproduces the measured values to high accuracy. 

In order to illustrate the impact of the <i-term correction in (}6"4*j) on the estimate of 
Yang-Mills scale parameters, we briefly consider the lattice scale parameter A[ at . This 
parameter is implicitly defined at the 2-loop level by considering (SO] 



a 



lat 



4ir I d' A 2 



lat 



In 



In 



« 2 A L 



(66) 



where for the SU (N) gauge group 



a 



lat 



2 -P. 

N H 



(67) 



Inverting eq. 



consistently up to 2-loop perturbation theory yields 



47T 



^ AL) = - f- 



a 



-i 

lat 



2fr /47T _! 



(68) 
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O gluon (no tadpole improvement) 




4 



- ln( a a ) - In (a c) 



Figure 2: The gluon and ghost renormalization constants, Z3 (/?,//) and Z 3 (/?,/i), for 
H = 3 GeV (left panel) y axis is arbitrarily scaled. Towards the cutoff dependence of 
Z\ (right panel). Figures corresponding to the Tubingen data only. 



Then, using a 



-1 

lat 



7T 



j3 from eq. ()67|) and eliminating a by subtracting eq. (|64|) from 



the latter equation we find 



'A 2 
In ' Alat 



lim 

/3->oo 



C — 



Thus, if we extrapolate to the continuum limit /3 



00 we obtain 



A 



lat 



e~ c/2 ^a = 0.112(5)^ 



(69) 



(70) 



Using the value a = [440 MeV] 2 one gets A lat = 49(2) MeV. If one instead of the limit 
in (|59*jl assumes that the asymptotic scaling regime is reached for, e.g., (3 = 2.5, one 
gets Ai a t = 0.0188(8) y/a. This is the order of magnitude familiar from the literature. 
Hence, for a scaling analysis of lattice results employing (3 e [2.15,2.8] the irrelevant 
term of order 1/(3 is still important. 



4 Results 

4.1 Renormalization constants Z3, Z3 and Z\ 

Let us firstly focus on the renormalization constants Z3 and Z%. As outlined in Subsec- 
tion l3.2l these constants are obtained from "matching" the lattice data from simulations 
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Figure 3: The analogue of Fig. |21 for the Sao Carlos data. The fits shown on the 
left-hand side for the gluon (upper curve) and ghost (lower curve) renormalization 
constants are done with 7 as a free parameter. All fits neglect the leftmost data point. 



using different ft values. Figures |2 and El show the cutoff dependence of these constants 
respectively for the Tubingen and Sao Carlos sets of data. 

Using the results above, we can check that our data are consistent with the pre- 
dictions from perturbation theory. For large enough UV cutoff, one expects that the 
1-loop behavior is recovered, i.e., 



Z a (J3,(i), Z 3 (j3,fi) « b [- \n(aa 2 )+u\ 



(71) 



If only small lattice spacings are considered, u is related to the ratio between the string 
tension and the Yang-Mills scale parameter at 1-loop level 



00 



In 



7r 2 a 



(72) 



loop 



Here, we treat u as a fit parameter and explore a range of lattice spacings where one 
would already expect significant deviations from the 1-loop behavior. As shown in 
Figures El and 01 (left panel), a good consistency with the known anomalous dimensions 
is observed. 

For the Tubingen data, we find that 10 w 1.13 is a good choice for reproducing the 
data for Z 3 and Z 3 simultaneously. It turns out that within the ft range explored in the 
Tubingen runs, tadpole improvement has a minor effect on the anomalous dimension. 
The adjoint plaquette used for the tadpole improvement [see eq. lJ33jl ] employing the 
adjoint representation is listed in Table 01 
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P 2.15 


2.2 


2.3 


2.375 


2.45 


2.525 


u a Q d P 0.225(42) 


0.241(43) 


0.274(45) 


0.297(46) 


0.317(46) 


0.336(47) 



Table 4: Expectation value of the adjoint plaquette (J2EJ used for the tadpole improve- 
ment of the gluon fields derived from the adjoint representation. 

For the Sao Carlos data we have performed the fits with 7 as a free parameter, 
leaving out the data point with (3 = 2.2 (i.e. the leftmost point in Figure EJ). We obtain 
the following values for the gluon and ghost cases 

7 g iuon = 0.60(5) and 7 ghost = 0.32(7) . (73) 

We see that the values are respectively consistent with 13/22 « 0.59 and 9/44 rs 
0.20 within error bars (but notice that there is a discrepancy of almost two standard 
deviations for 7 g host)- In this case we have not succeeded in finding a value of uj 
describing the behaviors for Z 3 and Z 3 simultaneously. 

In order to interpret the product of ghost and gluon form factors as the running 
coupling strength (see Subsection 12. 7|) . it is of great importance that the ghost-ghost- 
gluon- vertex renormalization constant Z\ be finite in the continuum limit. For detect- 
ing the UV behavior of Zx, let us investigate the product 

ZsG9, Ai) Z 3 2 (/3, /x) = z*(p, , (74) 

ao{A uv ) 

where (JSHJ) was used. The left-hand side of the latter equation can be directly obtained 
from the numerical result for the renormalization constants Z^(/3, //) and Z^(f3, fi). Note 
that for large UV cutoff, one finds 

1 A 2 

— — r oc In , = uv = — ln(aa 2 ) + constant, (75) 

ao(Auv) At loop 

where the constant comprises cutoff (and therefore 0) independent terms. The crucial 
point is that if the product Z 3 (fi, fi) Z%([3, /i) rises linearly with — ln(cra 2 ) the additional 
factor Zf((3,fi) must be finite in the continuum limit [since the renormalized coupling 
«ij(/i 2 ) is assumed finite]. Our numerical findings for Z 3 (P,/i) Z\{j5,\i) are also shown 
in Figures El and El (right panel). The data nicely support Taylor's findings, i.e. Z\ is 
cutoff- and therefore /i-independent. 

4.2 The running coupling constant 

Once it is established that Z\ is finite, the momentum dependence of the running 
coupling constant can be simply derived from the product (JBTj) 

a R (p 2 ) = a R (fi 2 ) F R (p 2 ,n 2 ) J 2 R (p 2 ,fi 2 ) . (76) 
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2-loop: A.„„ = 0.95GeV, p =2 GeV 

^ MOM ' 

I -loop contribution 



a 




123456 02468 10 12 

P [GeV] p [GeV] 

Figure 4: The running coupling in comparison with the results from perturbation 
theory: Tubingen data (left panel) and Sao Carlos data (right panel). The momentum 
cutoff is pm = 2 GeV in the former case and pu = 2.5 GeV in the latter. 

The overall normalization factor can be obtained by comparing the lattice data with the 
well-known perturbative result, which is valid at high momentum. Here, we compare 
with the 2-loop expression, which is known to be independent of the renormalization 
prescription, i.e., 

( 2 /A 2 \ 47r Ji 2/3 1 ln(lnx)\ 

a 2 . loop (x = p /A 2 _ loop ) = _|1 - j , (77) 

with j3 and Pi given in eq. (JoTj) . In order to obtain A 2 _; oop and to fix the overall 
factor, we fitted a#(/i) F R (p 2 , /i 2 ) Jj,(p 2 ,/i 2 ) to the 2-loop running coupling a 2 _; oop (p) 
where only momenta p > pu were taken into account. Fitting parameters were «_r(/z) 
and the 2-loop perturbative scale A 2 _j oop . Starting from a very low value pu we fit 
these parameters while gradually increasing pu- For small values of pm, we observe 
that the functional form of (fTTj) tries to incorporate genuine non-perturbative effects 
by adjusting A 2 -i oop , thus, introducing a spurious pu dependence to A 2 -z 00 p- However, 
a plateau is reached for the Tubingen data at pu ~ 2 GeV indicating that the data are 
well reproduced by the 2-loop formula in this regime. We find in this case that 

A2-ioo P = 0.95(15) GeV . (78) 

For the Sao Carlos data we have cut the data at pu ~ 2.5 GeV, which corresponds to a 
large drop in the x 2 / d.o.f. of the fit. (This also corresponds to reaching a relatively good 
plateau for A 2 _i op obtained from the fit.) We obtain the value A 2 _; oop = 1.2(1) GeV. 
The two values are consistent within error bars. 

Our final results for the running coupling constant are presented in Figure |U In 
the IR region, the two data sets show a clear departure from the perturbative behavior 
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Figure 5: The gluon form factor F^(p) and the ghost form factor Jr{p) as a function 
of the momentum transfer p (Tubingen data). 

and suggest a finite value a c for the running coupling constant at zero momentum. 
We estimate a c = 5(1). This value is in agreement with our previous fits for aji(p 2 ) 
jSHEi], and is consistent with the DSE result of Ref. a c ~ 5.2, using k = 0.5 as 
input. 

4.3 Gluon and ghost form factors 

Figure 03 shows the Tubingen data for the gluon and the ghost form factors. The error 
bars comprise statistical errors only. It turns out that one observes an additional scat- 
tering of the data points which is not of statistical origin. This additional systematic 
noise is pronounced when simulated annealing is used for gauge fixing and it afflicts 
especially the small momentum range. We attribute this error to the residual uncer- 
tainty of gauge fixing (Gribov noise). In particular, since the simulated annealing is 
capable of hopping from one local minimum to the other, the algorithm is sensitive 
to the large-scale structure of the minimizing functional. This hopping then produces 
a non-Gaussian noise which is underestimated when one uses the standard Gaussian 
error propagation. For this reason, we dropped the first three momentum points from 
the Tubingen data sets. 

In Figs. I^and [7| we report the rescaled Sao Carlos data for the renormalized gluon 
(respectively ghost) form factor and the data for the corresponding propagators. We 
stress that in the gluon case finite-size effects depend on whether we consider the full 
propagator or the form factor. In fact, for the propagator these effects are larger in 
the IR region, while for the form factor the effects are larger in the UV limit (due to 
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Figure 6: The gluon form factor Fr(p) (above) and the gluon propagator Dr(j>) (below) 
as a function of the momentum transfer p (Sao Carlos data). Note the logarithmic scale 
in the second plot. Fit of gluon propagator using the analogue of eq. (J75j) has been 
done for momenta p > 2 GeV and with A = 1.2 GeV. 



24 



Q. 



1 1 1 1 


1 1 






- + 


- 






- t 




\ 




+-H-HH-H- ++ ^_ 

1 1 1 1 


t- + + + ++ 

1 1 



2 4 6 8 10 12 

p [GeV] 




Figure 7: The ghost form factor Jr(j>) (above) and the ghost propagator Gr(p) (below) 
as a function of the momentum transfer p (Sao Carlos data). Note the logarithmic scale 
in the second plot. Fit of ghost propagator using the analogue of eqs. (fTTJj) has been 
done for momenta p > 2 GeV and with A = 1.2 GeV. 
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the multiplication by p 2 ). Thus, the ranges of momenta (for each 0) considered for the 
plots are different in the two cases. Nevertheless, the matching factors obtained are in 
agreement. The difference in finite-size effects between propagator and form factor is 
less pronounced when considering the ghost propagator. 



At sufficiently high momentum p > Pm (we found in Subsection 14.21 that pm ~ 
2 GeV is an acceptable choice), the momentum dependence of the renormalized form 
factors should be given by the formula 



F R (p 2 ,fi 2 



■/*(pV) 



d 2 {n) 



&2-loop 



A 2 

lv 2- 



loop 



1 + ^a 2 ~ioop(y 



A 2 



loop 



(79) 



where 7 is the leading-order anomalous dimension of the gluon (respectively ghost) 
propagator, given by 7 = 13/22 (respectively 7 = 9/44). The parameter 7 stems from 
the next to leading order to the anomalous dimension and is scheme-dependent. At 
least in the MS scheme, this parameter is small. Furthermore, Equation ()79)1 can 
be derived from the renormalization-group equation using the 2-loop scaling functions 
(3(g r ) and 7a(5 , b)- Hence, (I79J) originates from the resummation of an infinite set of 
2-loop diagrams and, e.g., comprises the so-called "leading logs". 

Using the Tubingen data set, pu = 2 GeV and A 2 -i 00 p = 950 MeV, we fitted 7 to 
the gluon and ghost data, respectively. We find that these parameters are indeed small, 
i.e., 

7giuon = -0.036(18) , 7 ghost = 0.011(10) . (80) 

Although the errors on these parameters are rather large, we find it encouraging that 
the parameters appear with opposite signs. In the case that the product of form factors 
Fr(p 2 ) Jji{p 2 ) is indeed renormalization group invariant, one would expect that 



7gluon ~\~ 2 7ghi 



ost 



0. 



It is clear from the Figures El and 03 that the high momentum tail is well reproduced 
by (EU). 

For the Sao Carlos data we have found that a cut at pu = 2 GeV corresponds to a 
large drop in the \ 2 /d.o.f. of the fits, both for the gluon and for the ghost cases. We 
start by fitting the leading-order term only (i.e. ignoring 7)- We get 



A 



2— loop 



1.19(4) and A 2 _ loop = 1.13(2) 



(82) 



respectively from the fits of the gluon and of the ghost propagator. These values are 
consistent with the result A 2 _i oop = 1.2(1), obtained in the previous section. We then 
fix this value for A 2 _i oop and perform fits with 7 as a free parameter. We obtain a good 
description of the data (see Figs. |B1 and HJ), with values for 7 even consistent with zero. 
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As can be seen from our plots, the gluon form factor is suppressed in the low 
momentum regime, while the ghost form factor is divergent. Correspondingly, the 
ghost propagator diverges faster than 1/p 2 and the gluon propagator appears to be 
finite. As mentioned in the Introduction, an IR- finite gluon propagator [2Ej [2] arid 
an IR-divergent ghost form factor [301 EE] were obtained before by separate studies. 

A quantitative analysis of the IR behavior for the propagators — including the 
evaluation of the exponent k mentioned in the Introduction — was already presented 
in Refs. jSHEHj- More thorough such analyses will be presented separately for the two 
sets of data in Refs. IMj. 

5 Conclusions 

For the first time, evidence from extensive lattice simulations is provided that the 
ghost-ghost-gluon-vertex renormalization constant Z\ is indeed finite in continuum 
field theory (as found by Taylor using all orders perturbation theory). Also, our result 
is probably not affected by Gribov ambiguities, since Z\ is obtained using data in the 
UV limit. It therefore appears that the Gribov ambiguities (in the lattice approach 
and the Faddeev-Popov Quantization) do not afflict the renormalization of the vertex. 

Also, we performed a thorough study of the gluon and the ghost form factors. Our 
data favor the scenario of an IR finite (or even vanishing) gluon propagator while the 
ghost form factor is singular in the IR limit. 

Finally, we have obtained the running coupling constant over a wide range of mo- 
menta using the data for gluon and ghost form factors. Our data are consistent with 
the existence of an IR fixed point a c = 5(1). Note that this value is inside the interval 
given by the DSE expression 

We stress that we compared our results for two slightly different lattice formulations, 
obtaining consistent results in all cases considered. 
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